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Abstract. Supersonic turbulence fragments the interstellar medium into dense sheets, 
filaments, cores and large low density voids, thanks to a complex network of highly 
radiative shocks. The turbulence is driven on large scales, predominantly by supernovae. 
While on scales of the order of the disk thickness the magnetic energy is in approximate 
equipartition with the kinetic energy of the turbulence, on scales of a few pc the 
turbulent kinetic energy significantly exceeds the magnetic energy. 

The scaling properties of supersonic turbulence are well described by a new analyt- 
ical theory, which allows to predict the structure functions of the density and velocity 
distributions in star-forming clouds up to very high order. 

The distribution of core masses depends primarily on the power spectrum of the 
turbulent flow, and on the jump conditions for isothermal shocks in a magnetized gas. 
For the predicted velocity power spectrum index f3 = 1.74, consistent with results of 
numerical experiments of supersonic turbulence as well as with Larson's velocity-size 
relation, one obtains by scaling arguments a power law mass distribution of dense cores 
with a slope equal to 3/(4 — /3) = 1.33, consistent with the slope of the Salpeter stellar 
initial mass function (IMF). Results from numerical simulations confirm this scaling. 
Both the analytical model for the stellar IMF and its numerical estimate show that 
turbulent fragmentation can also explain the origin of brown dwarfs. The analytical 
predictions for the relative abundance of brown dwarfs are confirmed by the observa- 
tions. 

The main conclusion is that the stellar IMF directly reflects the mass distribution 
of prestellar cores, due predominantly to the process of turbulent fragmentation. 



1 Introduction 



Turbulence in the interstellar medium (ISM) of the Milky Way - and more 
generally turbulence in the discs of other galaxies - is of crucial importance for 
both the structure and evolution of the galaxy. The importance of turbulence is 
both direct, through its influence on the pressure equilibrium and stratification, 
and indirect, through its influence on the star formation process. 

The vertical pressure equilibrium and stratification of the ISM is determined 
by the level of turbulence, together with the temperature distribution of the 
medium (which is in turn probably tightly coupled to the turbulence), and it 
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is likely that even the distributions of magnetic fields and cosmic ray particles, 
which also contribute to the pressure and stratification, are integral parts of the 
same process; it is unlikely that the near equipartition of the energy content of 
turbulence, magnetic fields, and cosmic ray particles is a mere coincidence. 

It has long been realized that turbulence in the interstellar medium, in partic- 
ular in the cold, molecular cloud components, is highly supersonic [ ^ , ^3|j44| , |26| , p7t . 
More recently, it has been realized that the supersonic nature of the turbulence 
is a boon, rather than a nuisance, when trying to understand the properties of 



the ISM, the cold molecular clouds, and star formation |59| , |64 25 . It turns out 
that supersonic turbulence is in many respect similar to ordinary, subsonic tur- 
bulence, and that it thus has a number of generic, statistical properties. Much 
like ordinary turbulence, its decay time is of the order of the dynamical time, 
even in the MHD-case p0p2||5l| , [79| , |64[ . And much like ordinary turbulence, it is 
characterized by power law velocity power spectra and structure functions over 
an inertial range of scales P, |TT| , |r2"] . 

An important difference between supersonic and subsonic turbulence is the 
distribution of density. A supersonic medium is, by definition, highly compress- 
ible; on average its gas pressure P g is small relative to the dynamic pressure pu 2 . 
As a consequence, a supersonic medium is characterized by a wide distribution of 
densities. A turbulent and isothermal supersonic medium has a log-normal den- 
sity probability distribution (PDF) j8|j6l0, with a dispersion of linear density 



proportional to the Mach number 




Cold molecular clouds arc indeed 



approximately isothermal, and are known to have a very intermittent density 
distribution, consistent with the properties of isothermal supersonic turbulence 
[ p0| . Deviations from isothermal conditions are in general of the type where com- 
pression leads to even lower temperatures (effective gas gamma less than unity) 
[ fri| , resulting in a density PDF skewed towards greater probability at high den- 
sities. The PDF may be described as a skewed log-normal, with a high density 
asymptote that formally tends to a power law in the limit T — > [f74]j5^] . 

Effectively then, supersonic turbulence acts to fragment the ISM, causing 
local density enhancements also over a range of geometrical scales. Molecular 
clouds themselves represent relatively large scale density enhancements, prob- 
ably caused by the random convergence of large scale ISM velocity features 
[pp|j33|] . Inside molecular clouds smaller scale turbulence leads to high contrast 
local density enhancements in corrugated shocks, intersections of shocks, and in 
knots at the intersection of filaments. Such small scale density enhancements are 
'up to grabs' by gravity; if their density is sufficiently high, relative to their tem- 
perature and the local magnetic field strength, they form pre-stellar cores, and 
eventually collapse to form stars. The decisive importance of turbulence in this 
process makes it possible to predict the distribution of masses of the pre-stellar 
cores, and hence the distribution of new borne stars, the initial mass function 
(IMF) HH). 

The process of star formation is indeed crucial to understand. Only by un- 
derstanding star formation, qualitatively and quantitatively, can we understand 
galaxy formation. We need to understand evolution effects to answer questions 
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such as "Was star formation different in the Early Universe?" . We need to un- 
derstand environmental effects to answer questions such as "Do other galaxies 
have different 'Larson laws' ?" 

We also need to understand star formation to answer questions related to 
Gamma- Ray Bursts; e.g., "Are Very Massive Stars progenitors of Gamma- Ray 
Bursts?" , and "What environment does the blast wave associated with Gamma- 
Ray Bursts encounter"? 

Finally, we would like to understand star formation as such, because it is 
a neat problem - one that involves supersonic, selfgravitating MHD turbulence 
and thus was thought to be enormously difficult. With access to supercomputer 
modeling the problem has become tractable, and it has turned out a posteriori 
that it is even partly amenable to analytical theory. 

In the subsequent sections of this tutorial star formation and turbulence in 
the interstellar medium is discussed in more detail. Section |^ discusses super- 
nova driving of the ISM, Section |^ discusses properties of supersonic turbulence, 
Section ^ summarizes a new theory of supersonic turbulence, while Section dis- 
cusses star formation and the initial mass function. Conclusions are summarized 
in Section |^. 

2 Supernova driving of the Interstellar Medium 

With turbulence being of such fundamental importance in determining the struc- 
ture and star formation efficiency of the interstellar medium it is important to 
understand what its primary sources are, and what its overall energy budget is. 

First, an estimate and lower limit of the energy input needed to sustain 
interstellar turbulence is given by Kolmogorov's scaling expression for the energy 
transfer rate in a turbulent cascade |39| , 

e ~ P U 3 /L , (1) 

where U and L are velocity amplitudes and length scales, respectively. In Kol- 
mogorov's classical theory this quantity is assumed to be invariant across the 
inertial range, and for our purposes this is adequate; subsequent enhancements of 
Kolmogorov's theory j75) and modifications for supersonic conditions Jl^] would 
not change the following estimates significantly. 

Observationally, the velocity dispersion in the ISM adheres to Larson's scaling 
law, 

U~ 1 kms- 1 (J^ , (2) 

with a w 0.4 p2p^ , p7[ |, which means that an estimate based on (^) only depends 
very weakly on the scale L on which the estimate is based. On scales L ~ I kpc, 
the turbulent velocity dispersion is of the order U ~ 10 kms -1 which 
leads to the estimate e ~ 5 10 50 ergkpc~ 3 Myr, using an average ISM density 
- 1.5 10~ 24 gem" 3 ||. 

For comparison, the rate of energy input to the ISM from supernovae is of 
the order of 10 53 ergkpc _3 Myr, based on a rate of one SN per 70 years in a 
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galactic volume spanned by a radius of 15 kpc and a disk thickness of 200 pc 
Thus, less than one percent of the average supernova energy input is necessary 
to sustain the turbulent cascade of energy in the ISM. 

Two questions come to mind: 1) Is there at all a turbulent cascade, and 2) 
is the energy from supernova at all available for feeding such a cascade? 

The answer to the first question is definitely affirmative; whatever its source, 
the observed velocity field at large scales can do nothing but drive a cascade 
towards smaller scale, since there is no dissipation mechanism that operates at 
such large scales. As has recently been shown |ll],[l2|], it makes little difference 
whether turbulence is subsonic or supersonic; similar cascades arise in both cases, 
only details such as power law exponents differ. The predicted scaling of the 
velocity dispersion with size is consistent with the observed (Larson's law) scaling 
[ ^2p^,p7t - Most of the observed scatter around the expected scaling (e.g., Fig. 
9 in |27[) is probably due to cloud-to-cloud variations - observations for a single 
cloud (Polaris) define a remarkably well defined scaling over more than three 
orders of magnitude in size p7| . A complementary piece of evidence for power 
law behavior comes from the observed relation between age difference and spatial 
separation |p5| . 

The answer to the second question is less obvious, but in the end also affir- 
mative. One might think that supernova energy input occurs at small scales, and 
hence cannot be a source at large scales for the turbulent cascade. However, as 
has been demonstrated by detailed numerical simulations fj^j32p f| , supernovae 
are indeed capable of sustaining a turbulent cascade with velocity dispersions 
consistent with observed values. The transfer of energy to large scales occurs 
through the expansion of supernova bubbles and super-bubbles; i.e., via the hot 
component of the ISM. The hot component coexists with cooler components (or 
rather a continuous distribution of temperatures) , and the expansion of the hot 
component into channels and chimneys creates kinetic energy on large scales, 
available for cascading to smaller scales, also in the cooler components. 

The numerical models demonstrate that supernova-driving of the interstellar 
medium is a viable and probably dominating mechanism, at least in the disc 
of our galaxy. The detailed numerical models are also broadly consistent with 
analytical and semi-analytical models of supernova feedback and turbulent self- 
regulation in galactic discs [ p0||8^ |. 

One may ask whether other sources of energy input could be significant. 
Winds from hot stars is one candidate that may contribute |85||39| ; regions that 
create supernovae of type II are likely to also contain hot, early type stars. Jets 
from new-borne stars have been mentioned as an energy input candidate, but it 
is unlikely to be significant on scales above a few pc. 

Larson's original paper lists irregularities and asymmetries in the rotation 
curve as larger scales that fit into a general power law behavior, at L ~ 2 — 10 
kpc, although with a break at L ~ 1 kpc. This could be taken as an indication 
that such irregularities, stemming for example from large scale density waves, 
could also be a source of driving for the turbulent cascade. But the relation could 
also go the other way; irregularities on scales of several kpc could be the imprint 



Star Formation and the Initial Mass Function 



5 



of old super-bubbles, stretched in the direction of rotation by the differential 
rotation. 

In other contexts, such as star-burst galaxies, the balance between contri- 
butions to the driving may be different; kinetic energy input from collisions or 
close interactions between galaxies may be an important energy input channel 
there, for example. 

3 Turbulent cascade of the Interstellar Medium 

When regarded as an isolated phenomenon, molecular clouds have traditionally 
given rise to concerns about the source of their turbulence, their life times, and 
about their support against gravity | ]76| , p3| . 

3.1 Molecular clouds as part of a turbulent cascade 

Molecular clouds are known to be significantly supersonic, with observed turbu- 
lent velocities of the order kms -1 , while typical sound speeds at molecular cloud 
temperatures are ~ 0.2-0.3 kms -1 . The supersonic velocities were assumed to 
give rise to very rapid dissipation in shocks, and hence explaining how the ob- 
served velocities are sustained was regarded as a problem. A popular suggestion 
for a solution to the problem was that the observed velocities are essentially 
magneto-hydrodynamic waves, with very low dissipation rates f|,|8^]. 

But in light of the conclusions of the previous section there is really no reason 
to be concerned about how the turbulence of molecular clouds is sustained; 
supply of kinetic energy at molecular cloud size scales is available from the 
turbulent cascade; i.e., simply from larger scale motions. In fact, a molecular 
cloud is probably borne precisely because the larger scale velocity field happens 
to have a local maximum of convergence there f|,|]j33) . 

3.2 Supersonic turbulent cascades 

Figure |l] shows examples of power spectra of supersonic turbulence, for two 
values of the numerical resolution. The numerical simulations were performed 
with the same method as in pl| , |6^ , |67]j6^ ] ; a fifth order in space, third order 
in time staggered mesh method, using per-unit-mass variables. The turbulence 
is driven by a random external force that is applied in 1 < k < 2 in Fourier 
space. Only the solenoidal components of the force are used. In order to ensure 
a smooth forcing the time derivatives of the Fourier components of the force 
are regenerated at time intervals of about one dynamical time and the force is 
computed from a time integral. 

An inertial (power law scaling) range is present, and extends to progres- 
sively higher wavenumbers at higher numerical resolution. How is it possible 
that supersonic turbulence gives rise to a turbulent cascade much like that of in- 
compressible turbulence? One clue comes from the dissipation rate of supersonic 
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Fig. 1. The solenoidal power spectrum, compensated by fc 1,74 , and the ratio of com- 
pressional to solenoidal power in M ~ 10 numerical experiments with resolution 250 3 
(diamonds) and 500 3 (stars), using random solenoidal driving at 1 < k < 2. Dashed 
lines show comparison slopes with spectral indices -1.64, -1.74, and -1.84, respectively. 



turbulence. Numerical experiments revealed it to be similar to that of incom- 
pressible turbulence, if expressed in terms of dynamical times Tdyn = 
even for MHD-turbulence ||,|5(|[52 51 6^,^] • In qualitative terms one reason why 
the dissipation rate is not as high as was naively expected is that shocks in three- 
dimensional supersonic turbulence are typically oblique rather than head-on, and 
that fragmentation decreases the efficiency of interaction |7^,^| . 

The ratio of compressional to solenoidal kinetic energy is small in isotropic 
supersonic turbulence; typically E c /E s ~ 0.1 — 0.2 [|ll) (cf. Fig. |]). To appreciate 
why this is so it is helpful to consider the velocities on either side of a shock sheet 
formed by two interacting large scale streams. By definition, the gas upstream 
of the two stand-off shocks on either side of the sheet have no casual connection, 
and their orientations are therefore random with respect to each other. 

In terms of a coordinate system with one axis perpendicular to the plane of 
the sheet, two of the three velocity components are parallel to the sheet and 
hence incompressible, and only one is perpendicular (compressional). The com- 
pressional component is the one that gives rise to the shock, with its associated 
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Fig. 2. Density in an arbitrary cut through a numerical experiment with M w 10 
turbulence at a resolution of 500x500x500. a) Linear scaling (printed with gamma=0.7), 
normalized to the local maximum of the density, b) Logarithmic scaling, normalized 
to the local maximum and minimum of the density. 

stagnant region inside the shock sheet. It follows that one should indeed expect 
the compressional component to, on the average, carry less than one third of the 
kinetic energy, consistent with what is found in the numerical experiments. 

More generally, one may think of three-dimensional supersonic turbulence as 
an ensemble of shock sheets, and their associated intersections in filaments and 
knots. The typical history of a trace particle in such a flow is that it participates 
in a series of oblique shocks where, in each shock, the particle looses some of its 
kinetic energy. 

If the system consisted of an ensemble of stationary, plane-parallel shock 
sheets a fluid parcel would first hit one sheet, where its perpendicular kinetic 
energy would be essentially lost. It would then slide along the sheet until it 
hit the filamentary intersection with another sheet, and finally slide along the 
filament until it ended up in the stagnant region of a knot-like intersection of 
filaments. 

In a more general picture shock sheets are neither stationary nor plane- 
parallel, which allows trace particles to participate in a more extended series of 
shocks. As its kinetic energy is gradually reduced, so is the scale over which new 
shocks are likely to be produced. 

Figure || illustrates the structure of the density field in supersonic turbulence, 
modeled at a resolution of 500 3 . The left hand side panel shows linear density. 
Due to the large density contrast, only a few shock sheets and filaments are 
visible. The right hand side panel shows logarithmic density, and illustrates the 
general presence of intermittent density structure over a range of scales and 
density levels. 

The history of a fluid parcel in the real ISM might be as follows: It achieves 
initial, large kinetic energy by either being part of the ejecta from a supernova 
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Fig. 3. Histograms of extinction (left panels) and plots of dispersion of extinction in 
cells versus the mean cell extinction (right panels) . The top panels show the result from 
a super-Alfvenic model while the middle panels are from an equipartition model Q]. 
The bottom panels are observational data for the cloud IC5146 [|fl| 
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Fig. 4. Histograms of integrated antenna temperature of synthetic CO spectra j64], 
for a super-Alfvenic model (thick line), for an equipartition model (thin line), and for 
the Perseus cloud (dotted line) 



or, more likely, by being hit by the ejecta from a supernova. It becomes further 
compressed as the stream to which it belongs collides with other streams. As 
density increases cooling becomes significant, and the temperature decreases. 
The parcel may eventually end up as part of a molecular cloud. Inside the cloud, 
the process repeats itself, on successively smaller scales, creating in the end a 
shock core massive enough to form a star. More likely, though, the fluid parcel 
ends up in a structure too small to collapse by self-gravity, where it survives 
until being hit by the blast wave from another supernova, or the wind from a 
new-borne massive star. 



3.3 Super-Alfvenic conditions 

The proposal that magneto-hydrodynamic waves are main contributors to the 
velocity field in molecular clouds is now obsolete for several independent reasons. 
First, with the velocity field of molecular clouds part of a turbulent cascade there 
is no longer a problem to sustain the motions. Second, with the demonstrations 
that MHD-turbulence decays more or less as rapidly as hydrodynamic turbulence 
[jsl |50| , |5^ , [51 64, 7§| , the presumed 'advantage' of MHD-turbulence has gone away. 



Third, with observational and theoretical evidence that star formation takes 
place on a time scale not much longer than a crossing time [^Bf the required life 
times of molecular clouds are much shorter than was assumed in earlier work. 

The notion that the velocity field in molecular clouds is essentially Alfven 
waves lead to the assumption that the kinetic and magnetic energy in molecular 
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Fig. 5. Scatter plots of equivalent width versus velocity integrated antenna tempera- 
ture for two observed regions (two upper panels), for a super- Alfvenic model, and for 
an equipartition mode. The diamond symbols show the mean value of the equivalent 
width in each interval of integrated antenna temperature and the "error bars" show 
the one a distribution around the mean 



clouds are in near equipartition. Although inverting observations to obtain the 
magnetic field strength, density and velocity in the same structures is notori- 
ously difficult (cf. the discussion in Section 4.1 of p0|), equipartition remained 
a popular null-hypothesis. 

With access to numerical simulations it is possible to use the more robust 'for- 
ward analysis' method, where synthetic diagnostics computed from the results 
of numerical simulations are compared directly with the corresponding obser- 
vational diagnostics. Comparisons of extinction statistics, synthetic molecular 
lines, the antenna temperature - line width relation, and the statistical upper 
envelope relation between density and magnetic field strength (cf. Figs. all 
lead to the same conclusion; models with equipartition between magnetic and 
kinetic energy are inconsistent with the observations while models where the 
kinetic energy dominates over the magnetic energy (super- Alfvenic models) are 
consistent with the observations [ |60||6^ , |63| , |6^| . 

A direct illustration of the consistency of super- Alfvenic conditions with Zee- 
man observations of magnetic field strength is given in Fig. [?]. Note that several 
cores with B in excess of 100 [iG are found, even though the average B in the 
simulation is only 2.4 (J.G. This is a good example of the power of forward com- 
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Fig. 6. The B — n relation: observations and theoretical models. The thick contour 
lines are from the a super-Alfvenic model and the thin contour lines are from an 
equipartition model [B4| 




Fig. 7. Magnetic field strength versus H2 column density. Asterisks represent Zeeman 
splitting measurements jl6| ] (detections and uppper limits) with a least squares fit 
(dashed line). Squares represent cores from a super- Alfvenic numerical experiment 
(average B = 2.4 \xG) with a least squares fit (solid line). The dotted-dashed line 
marks equality between magnetic and gravitational energies 
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Fig. 8. Non-thermal line width of observed (squares, with dashed line least squares fit) 
and modeled (stars, with solid line least squares fit) NH3 cores versus their size |37 
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Fig. 9. Rotational velocity (upper panel) and velocity dispersion (lower panel) for 
numerical cores. 
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parisons in situations with strong intcrmittency; it would be very difficult to 
recover the mean field strength, or the mean magnetic energy, directly from the 
observations, which sample only the very small fraction of the cloud volume filled 
by the densest regions. Further illustrations are given in Figs. which show 
comparisons of velocity statistics with observations. 

Figure ^| is a comparison of the correlation of non-thermal line width with 
size in NH3 cores from the compilation by Jijina, Myers & Adams |3?J and in 
cores selected from a simulation of supersonic and super-Alfvenic turbulence ]67| . 
The least squares fit to the observational data yields the power law exponent 
0.56 ± 0.22 and the fit to the data the exponent 0.57 ± 0.15. Figure ^ shows the 
correlation of rotational velocity (upper panel) and internal velocity dispersion 
(lower panel) with size, for the numerical cores only. The rotational velocities 
are very low and of the order of the sound speed, as found in the observational 
data. In both panels of Fig. |^ the asterisk symbols that correspond to a size 
of almost 3 pc provide the values of rotational velocity and velocity dispersion 
computed over the whole simulated volume. Although the least squares fits are 
computed only for the cores, the values for the whole system are consistent with 
the fits. 



3.4 The magnetic flux problem 

A "magnetic flux problem" is often mentioned in this context. It is argued that 
since the average density of molecular clouds is at least 100 times larger than 
the average density in the galactic disk, and assuming magnetic flux freezing, 
the average magnetic field strength of molecular clouds should be much larger 
than the average galactic values of a few p.G. As demonstrated by Fig. [?] there 
is in fact no real problem - observations of core magnetic fields are completely 
consistent with predictions from models with average magnetic field strengths 
of a few \xG - there is at most a conceptual / perceived problem. 

This conceptual problem has a straightforward solution, which, ironically, is 
most easily demonstrated by the equipartition model. It has been shown in many 
numerical works that supersonic turbulence, even with equipartition of kinetic 
and magnetic energy (the traditional model for molecular clouds) generates a 
complex density field, with very large contrast sheetlike and filamentary density 
structures. These density enhancements do not correspond to significant varia- 
tions of the magnetic field strength in equipartition models, since in them strong 
compression can occur only along magnetic field lines. To the extent that tur- 
bulence on large scales (disk thickness) has approximate equipartition of kinetic 
and magnetic energy, molecular clouds can still easily form, as a consequence of 
compressions along magnetic field lines. 

Once a cloud is formed by large scale equipartition turbulence, it has a mean 
magnetic field strength close to the galactic value, and its internal dynamics 
is super-Alfvenic, because of the much increased density. Equipartition on the 
large scale, therefore, is not a problem for the origin of super-Alfvenic clouds. 

The argument applies recursively; should the super-Alfvenic cloud by chance 
create a region with local equipartition, further increase of the density is still 
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possible, by inflow of mass along magnetic field lines. One sees the statistical 
outcome of this in the B-n relation (Fig. ||); for any given density there is a 
wide distribution of magnetic field strengths, up to an upper envelope given by 
approximate equipartition. 

Inflow along magnetic field lines is also likely to occur in the phase when 
gravitation has taken over, after local cores are formed along filaments and in 
corrugated sheets. In that situation the magnetic field has already been com- 
pressed, and is oriented predominantly along the same filaments and sheets that 
also contain abundant mass at high density. 



3.5 Gravitationally bound and unbound clouds 

The scenario where star formation takes place in essentially a crossing time [^5| 
also alleviates earlier concerns about how to support molecular clouds against 
gravitational collapse. The gravitational binding energy of molecular clouds is 
often comparable to their turbulent kinetic energy, and hence quite a bit larger 



than their thermal energy |43[|27|| . This raised the question of the support of the 
clouds against gravitational collapse. Could the turbulent velocities be translated 
into a turbulent pressure that was able to support the clouds against collapse 
[ fu||l4| , |84| , [38[ ? Or did the solution lie in the observed, strong fragmentation of 
the medium @|^? 

In the 'turbulent fragmentation & star formation in a crossing time' paradigm 
it is natural to find some molecular clouds with roughly virial mass, as well as 
some with substantially less than virial mass, while there are essentially none 
with much larger than virial mass (cf. Fig. 7 of ^7|). Turbulent fragmentation 
creates clouds, initially without regard to gravity. Some of the clouds that are 
produced are gravitationally unbound (but may contain sub-structures that are 
gravitationally bound). Some other clouds are massive enough to be gravita- 
tionally bound (at least until their first supernovae blow out a major part of the 
cloud gas). Clouds that are created with a mass larger than virial start to col- 
lapse, which increases the velocity dispersion of their sub-structures until they 
appear to be essentially virial. 

The latter case represents the most direct and simple mechanism by which 
turbulence prevents global collapse of molecular clouds; i.e., through fragmen- 
tation rather than through "turbulent pressure" . The mechanism may be illus- 
trated by considering that extreme intermittency caused by strongly supersonic 
turbulence and cooling could create conditions where individual density max- 
ima move in essentially ballistic orbits relative to one another (73). Even under 
less extreme conditions intermittency may cause individual density maxima to 



collapse, while the cloud as such does not |59|,|61|,|38p 



3.6 Power laws and equipartition 

As mentioned above, even highly supersonic turbulence is characterized by power 
laws [pd|Jr^|. However, because of the strong intermittency of density and its 
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correlation with the velocity field, the spectrum of kinetic energy is not the 
same as the power spectrum of velocity. 

It is appropriate to define the spectrum of kinetic energy as the power spec- 
trum of p 1 / 2 ^, since the sum of squares of its Fourier components is equal to 
the kinetic energy. Empirically, from numerical simulations, one finds that the 
spectrum of kinetic energy is quite a bit more shallow than the power spectrum 
of velocity. The latter has a power exponent consistent with the theoretical ex- 
pectation (3 = 1.74 The former has a power exponent /?k. c . ~ 1-1. 

The power spectrum of the magnetic field is approximately parallel to that 
of velocity in the high-k part of the inertial range, and hence the spectrum of 
magnetic energy is steeper than the spectrum of kinetic energy. This may appear 
strange, at first. Why would the magnetic field have a power spectrum similar 
to that of velocity, when magnetic energy, ^B 2 , is measured in the same units 
as kinetic energy ^pu 2 and not in the units of velocity power u 2 l A possible 
explanation is that B 2 is weighted more towards the kinetic energy of the bulk 
of the volume. Assuming a log- normal PDF of density, with a dispersion of linear 
density J55| 

o P = bM , (3) 
where b rs |, the most common density is 

(1 + b 2 M 2 ) 2 

which is smaller than the average density (p). Thus, if the magnetic energy is 
in equipartition with the kinetic energy at those densities, rather than at the 
higher densities towards which the average kinetic energy is weighted, this would 
explain both why the spectrum of B is similar to that of u and why the average 
magnetic energy is below equipartition with the kinetic energy. 

With a difference in the power law exponents a gap develops from the (ob- 
served) equipartition at large scales (~ 100 pc). A power law index difference of 
0.65 implies E mag / Eki n < 0.1 at pc scales, which is consistent with the forward 
analysis of numerical simulations [6^ ] . 

Note that the discussion above is complementary to the one at the end of 
Sect. |3.3| - both views are helpful for understanding why small scale ISM motions 
are super-Alfvenic. 



4 A new analytical theory of supersonic turbulence 

Due to the complexity of the Navier-Stokes equations, mathematical work on 
turbulence is often inspired by experimental and observational measurements. 
Since geophysical and laboratory flows are predominantly incompressible, tur- 
bulence studies have been limited almost entirely to incompressible flows (or to 
infinitely compressible ones, described by the Burgers equation). Little attention 
has been paid in the past to highly compressible, or super-sonic turbulence. 
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Turbulent flows are traditionally described statistically by the structure func- 
tions of their velocity field (28) . The structure functions are defined as 



where u is the component of the velocity field perpendicular (transversal struc- 
ture functions) or parallel (longitudinal structure functions) to the vector t. 
In the inertial interval the structure functions obey scaling laws and the ex- 
ponent ((p) can be determined. The power spectrum of the velocity is the 
Fourier transform of the second order structure function, and may be expressed 



One may think that the study of high order structure functions is interesting 
only for testing models of intermittency in turbulent flows, and not very useful 
in the context of ISM turbulence and star formation. Actually, the intermittent 
nature of turbulence is crucial in modeling the process of star formation driven 
by turbulent fragmentation. Stars are formed in the densest regions of turbulent 
flows. These regions contain only a few percent of the total mass and fill an almost 
insignificant fraction of the total volume of a star forming cloud. High order 
moments defining the tails of statistical distributions of velocity and density 
are therefore very important in the process of star formation. Furthermore, low 
order density structure functions, which are obviously important to describe 
basic properties of turbulent fragmentation, can be shown to depend on velocity 
structure functions of very high order [p"2|| . 

The scaling of the velocity structure functions in incompressible turbulence 
is best described by the She-Leveque formula j75[ , 



The scaling exponents are computed relative to the third order, ((p)/((3), be- 
cause according to the concept of extended self-similarity |7],[l9| the relative ex- 
ponents are universal and better defined than the absolute ones. 

Boldyrev Q has proposed an extension of the She-Leveque's formalism |75| 
to the case of supersonic turbulence. Based on the physical interpretation of 
© by Dubrulle |l9[ ] , a fundamental parameter in the derivation of the velocity 
structure functions is the Hausdorff dimension of the support of the most singular 
dissipative structures in the turbulent flow. In incompressible turbulence the 
most dissipative structures are organized in filaments along coherent vortex tubes 
with Hausdorff dimension D — 1, while in supersonic turbulence dissipation 
occurs predominantly in sheet-like shocks, with Hausdorff dimension D = 2. The 
new velocity structure function scaling proposed by Boldyrev || for supersonic 
turbulence is 



S p (£) = (\u{x + £)- u(x))\ p ) cx 



(5) 



as E(k) oc k-f 3 oc fc- 1 -^ 2 ). 







This velocity scaling has been found to provide a very accurate prediction for 
numerical simulations of supersonic and super- Alfvenic turbulence 111]] , and has 
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been used to infer the structure of the density distribution in turbulent clouds 
& 

5 Star formation and the Initial Mass Function 

At least three unrelated ways of explaining the process of star formation and 
the origin of the stellar initial mass function (IMF) may be found in the litera- 
ture: i) Ambipolar drift contraction of sub-critical cores [[76p[|; ii) opacity-limited 
gravitational fragmentation p^gl^ll^gl^Q ; and iii) turbulent frag- 
mentation ign^mmggi . 



The first type of models rely on the assumption that both protostellar cores 
and their parent clouds are long lived systems in near equilibrium, supported 
against their gravitational collapse by magnetic field pressure. As discussed 
above, this assumption has been proven incorrect based on observational data 
and is inconsistent with the turbulent nature of star-forming clouds p^ , p5p^ j65[| . 
Furthermore, these type of models do not address the problem of the formation 
of massive stars or brown dwarfs, and have traditionally focused more on the 
evolution of individual protostars, without providing a self-consistent picture for 
the origin of the initial conditions. 

The second type of models is also inconsistent with the properties of star- 
forming clouds, because it applies the concept of gravitational fragmentation to 
the large scale, in the attempt of modeling the formation of a whole stellar popu- 
lation. The concept of gravitational instability is based on a comparison between 
the gas thermal and gravitational energies to define the smallest unstable mass, 
or Jeans' mass. However, star-forming clouds, as any region of the cold ISM 
above a scale of approximately 0.1 pc, contain a kinetic energy of turbulence 
that is much larger (typically 100 times larger) than their thermal energy, mak- 
ing the comparison of thermal and gravitational energies irrelevant on the large 
scale. Attempts to redefine the Jeans' mass Jl|,|| assuming that turbulence can 
provide pressure support against the gravitational collapse are flawed, because 
they miss the basic point that supersonic turbulence is actually fragmenting the 
gas. The main effect of the large kinetic energy of turbulence, relative to the ther- 
mal energy, is that the gas density and velocity fields in star-forming regions are 
highly non-linear, against the assumption of the gravitational instability model. 
In other words, clouds are already fragmented by turbulence, quite independent 
of their self-gravity. 

The third type of models, which we refer to as turbulent fragmentation mod- 
els, focus on the importance of the observed supersonic turbulence in molecular 
clouds and are therefore consistent with the large scale dynamics of star-forming 
regions. The idea of star formation driven by supersonic turbulence was pro- 
posed twenty years ago by Larson but has become popular only in the 
last few years, thanks to the progress of numerical simulations of supersonic 
magneto-hydrodynamic (MHD) turbulence. 

According to the model of turbulent fragmentation, protostellar cores are 
formed from gas compressed by shocks in the supersonic turbulent flow Ipq]. 
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While scale-free turbulence generates a power law mass distribution down to 
very small masses, only cores with a gravitational binding energy in excess of 
their magnetic and thermal energy can collapse. The shape of the stellar IMF is 
then a power law for large masses, since the majority of large cores are larger than 
their Jeans' mass. At smaller masses, the IMF flattens and then turns around 
according to the probability of small cores to be dense enough to collapse, which 
is determined by the PDF of gas density. 

5.1 The Initial Mass Function 

The mass distribution of dense cores formed in a supersonic turbulent flow follows 
from a number of properties of such flows J66| : i) The power spectrum of the 
turbulence is a power law; ii) the dynamics on scales covered by the power law 
is approximately sclfsimilar; iii) the typical size of a dense core scales as the 
thickness of the postshock gas; iv) the relevant shock jump conditions are those 
of MHD shocks. 

These properties are to a considerable extent already verified by numerical 
simulations, which also produce corresponding numerical IMFs $57\ , but for the 
purpose of deriving a theoretical IMF one may just adopt them as assumptions 

& 

The second assumption, about approximate selfsimilarity, is a crucial one. 
Since velocity amplitudes do depend on scale (by the first assumption and by 
Larson's relation) it is not an exact property. Nevertheless, since these flows 
are supersonic they follow essentially inertial paths in a large fraction of space 
(upstream of shocks). The thickness and density of the downstream, shocked gas 
does depend on the Mach number, and hence on the scale, but the filling factor 
of the shocked gas is quite small and does not disturb the overall selfsimilarity 
much. 

The distribution of cores that form in the shocked, downstream gas may be 
regarded as a distribution over linear sizes L of the upstream flows out of which 
they formed. By the assumption of approximate selfsimilarity the number of such 
regions per unit log L scales as L~ 3 . The upstream (Alfvenic) Mach number is 
denoted A4(L) and is assumed to scale as L a , where a is related to the power 
spectrum index /? = 1 + 2a. 

The typical mass of the cores that form in the shocked gas scales as X 3 pi, 
where (from the MHD shock jump conditions) A ~ Lj M.{V) is the thickness of 
the postshock gas, p\ ~ poA4(L) is its density, where po is the upstream mass 
density (similar to the mean density). The typical core mass is thus 

m(L)~ P0 L 3 /M(L) 2 ~L 3 - 2a , (8) 

which leads to the following expression for the mass distribution of dense cores: 

N(m) ~ L(mf ~ m" 3 /^ 2 ") ~ m - 3 '^ . (9) 

If the power spectral index (3 is consistent with the observed velocity dispersion- 
size Larson relation |^3| and with the numerical and analytical results P]lT||, 
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then (3 w 1.74 and the mass distribution is 

N(m) dlogm oc m _1 ' 33 dlogm , 



(10) 



which is almost identical to the Salpeter stellar IMF |72| ]. 

If L is defined as the scale of a molecular cloud, with average mass density 
po and Alfvenic Mach number Mo, the mass of the largest cores formed by 
turbulent fragmentation is estimated to be 



PoLq 



11 -max a j9 ' 

Ml 

10 4 M R and Mach number Mq ~ 10, m„ 



(11) 



In MCs with mass Mq 
100 M Q . 

While the majority of massive cores are larger than their Jeans' mass, mj, 
the probability that small cores are dense enough to collapse is determined by the 
PDF of the density of the cores, which is approximately Log-Normal. Even very 
small (sub-stellar) cores have a finite chance to be dense enough to collapse. 
If p(mj) dm j is the Jeans' mass distribution obtained from the PDF of gas 
density |3l| , the fraction of cores of mass m with gravitational energy in excess 
of their thermal energy is given by the integral of p(mj) from to m. The mass 
distribution of collapsing cores is therefore 



N(m) dlogm oc m" 3/(4 - /3) 



p(mj)dmj 



dlogm 



(12) 



The mass distribution is plotted in Fig. 10, for (3 — 1.8. In the top panel the 
mass distribution is shown for three different values of the largest turbulent scale 
Lq, assuming Larson type relations Ei| to rescale the average gas density, (n), 
and the rms Mach number, M, as a function of size, Lq. The mass distribution 
is a power law, determined by the power spectrum of turbulence, for masses 
larger than approximately 1 m Q . At smaller masses the mass distribution flat- 
tens, reaches a maximum at a fraction of a solar mass, and then decreases with 
decreasing stellar mass. The mass distribution peaks at approximately 0.4 m Q 
for the values M = 10, (n) = 500 cm" 3 , T = 10 K and /3 = 1.8, typical of 
nearby molecular clouds. Collapsing sub-stellar masses are found, thanks to the 
intermittent density distribution in the turbulent flow. This provides a natural 
explanation for the origin of brown dwarfs. 

Note that the power law shape of the IMF for mass values larger than about 1 
m Q is not affected by the average physical properties of the system. On the other 
hand the abundance of brown dwarfs is very sensitive to the average gas density 
and the rms Mach number of the flow. The middle and bottom panels of Fig. fLO] 
show the dependence of the mass distribution on the rms Mach number of the 
flow and on the average gas density respectively. One can see in the middle panel 
that for an average gas density of (n) = 500 cm~ 3 and an rms Mach number 
M = 5, typical of a molecular cloud complex such as Taurus, brown dwarfs 
are very rare, while for the same average gas density and an rms Mach number 
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Fig. 10. Mass distribution of gravitationally unstable cores from equation (|l2|). Top 
panel: Mass distribution as a function of the largest turbulent scale La, assuming 
Larson type relations (for rescaling (n) and M with Lo), To = 10 K and /3 = 1.8. 
Middle panel: Mass distribution as a function of the rms Mach number of the flow, 
assuming (n) = 500 cm -3 , To = 10 K and j3 = 1.8. Bottom panel: Mass distribution 
as a function of (n), assuming M = 10, To = 10 K and j3 = 1.8. 
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M. = 10, typical of a molecular cloud complex such as Orion (the density may 
be even larger), brown dwarfs are very abundant (even more abundant if the 
IMF were plotted in units of linear mass interval). This prediction is in fact 
unambiguously confirmed by the observations |47| . 

The thermal Jeans' mass is a more strict condition for collapse than the mag- 
netic critical mass. The magnetic critical mass depends on the core morphology 
in relation to the magnetic field geometry and strength. The latter correlates 
with the gas density with a very large scatter |34|. It is possible therefore that 
magnetic pressure support against the gravitational collapse limits the efficiency 
of star formation, while its effect on the shape of the mass distribution is of 
secondary importance. 

Observations show that the stellar IMF is a power law above 1-2 m 9 , with 
exponent around the Salpeter value x — 1.35, roughly independent of environ- 
ment p2|,p4|, gradually flattens at smaller masses, and peaks at approximately 
4|Q. The shape of the IMF below 1-2 m Q , and par- 
ticularly the relative abundance of brown dwarfs, may depend on the physical 
environment (47). These observational results are all consistent with our theo- 
retical IMF. 

It has been argued that only a small fraction of the mass of each collapsing 
core may end up into the final star, due to mass loss in protostellar winds, with 
a major effect on the stellar IMF. However, stellar winds could be important 
for the origin of the stellar IMF only if the ratio of initial core mass to final 
stellar mass were comparable to the total mass range for stars (~ 10 4 , from 
~ 100 M Q to ~ 0.01 Mq), as pointed out by Elmegreen |||. This is highly 
unlikely, because i) the correct slope and mass range of the IMF is already 
achieved by turbulent fragmentation alone and ii) observational results indicate 
that the mass distribution of prestellar cores is indistinguishable from the stellar 
IMF |5^,^,^2|,Q , as predicted in earlier work on turbulent fragmentation and 
the origin of the stellar IMF Q . 



0.2-0.6 m Q HP,p 46 



5.2 Mass distribution of prestellar cores in numerical simulations 

The mass distribution of prestellar cores may be measured directly in numerical 
simulations of supersonic turbulence. With a mesh of 250 3 computational cells, 
and assuming a size of the simulated region of a few pc, it is not possible to follow 
numerically the gravitational collapse of individual protostellar cores. However, 
dense cores at the verge of collapse can be selected in numerical simulations by 
an appropriate clumpfind algorithm. Such an algorithm should scan all density 
levels and recognize when a large core is fragmented into smaller and denser 
ones, in which case the large core should not be counted. Cores should also be 
excluded if their gravitational energy is not large enough to overcome thermal 
and magnetic support against the collapse, since only collapsing cores should be 
selected. 

A mass distribution of collapsing cores, derived from the density distribution 
in a numerical simulation is shown in Fig. O. The computational box with 128 3 
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Fig. 11. Mass distribution of collapsing cores in the range 0.2-100 Mq, derived from 
a 128 3 numerical simulation. 



cells has been scaled to two scale ranges, suitable for sampling cores in the inter- 
vals 0.2-2 M and 2-100 M , respectively. The mass distribution above 1 M Q 
is a power law consistent with our analytical result and with the observations. 
Below 1 Mq the histogram flattens and then turns around at approximately 
0.3 Mq, also consistent with the analytical theory and the observations. The 
cut-off at ~ 0.2 M Q is due to the finite numerical resolution; the grid size, rms 
Mach number, and mean density together impose a limitation on the mass of 
collapsing cores. 

Stretching the mass interval of sampled cores further into the brown dwarf 
regime requires larger numerical resolution. Figure [l^ shows the mass distribu- 
tion of collapsing cores derived from two snapshots of a 250 3 simulation. The 
average gas density has been scaled to 500 cm -3 and the size of the computa- 
tional box to 10 pc. These values have been chosen to be able to select cores 
in a range of masses from a sub-stellar mass to approximately 10 Mq. With 
this particular values of average gas density and size of the computational box, 
the smallest mass that can be achieved numerically is 0.057 Mq. Brown dwarfs 
masses (< 0.08 Mq) are therefore included. With an even larger numerical mesh, 
or assuming a larger average density (and a smaller size), even smaller masses 
would be selected. Turbulent fragmentation thus provides a natural explanation 
for the origin of brown dwarfs. This was found from the analytical model of the 
IMF presented above and is here confirmed from the numerical mass distribution. 
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Fig. 12. Solid line: Mass distribution of collapsing cores, derived from the density 
distribution of two snapshots of a 250 3 simulation with rms Mach number M ~ 10. 
The simulation is scaled to physical units assuming (n) = 500 cm -3 , To = 10 K, and a 
size of 10 pc. The fractional mass in collapsing cores is 5% of the total mass. Dashed 
line: Analytical mass distribution, N(m), computed for (n) = 500 cm -3 , To = 10 K, 
(3 = 1.8. 



Observed star-forming clouds appear very filamentary, as the projected den- 
sity field of supersonic turbulent flows. We have performed accurate comparisons 
of statistical properties of turbulent flows with observational data, by computing 



synthetic spectral maps of molecular transitions |62 63.65|] . The synthetic spec- 
tral maps are obtained by computing the non-LTE radiative transfer problem 
using the density and velocity fields of the MHD simulations. We have shown 
that fundamental statistical properties of supersonic turbulence are unambigu- 
ously found in the observational data of star-forming clouds |35| . 

In star-forming clouds, prestellar cores and young stars tend to concentrate in 
the densest filaments and cores. Since filaments and cores of the same nature are 
found in the numerical simulations, it is interesting to visualize the position of the 
collapsing cores selected numerically, relative to the gas density distribution. In 
Fig. [H^a voxel projection of the density field is shown, where all the numerically 
selected cores have been highlighted as bright spheres. The size and brightness 
is a function of the core mass. The brightness is also a function of the optical 
depth of the gas between the observer and the stars, to mimic the effect of dust 
extinction. Figure O shows beautiful filamentary structure in both the gas and 
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Fig. 13. Voxel projection of the density field of a snapshot of a 250 3 numerical sim- 
ulation of supersonic and super-Alfvenic turbulence. Collapsing cores are highlighted 
as bright spheres, with brightness and size varying as a function of the core mass. 
The brightness also depends on the column density of gas between each core and the 
observer, in order to mimic the effect of dust extinction. Slightly extended patches of 
bright emission are "unresolved" stellar clusters. The fraction of mass in collapsing 
cores is about 5% in this simulation. 

the stellar distribution, very reminiscent of observed star-forming regions. It is 
quite amazing that a numerical simulation of randomly driven supersonic and 
super-Alfvenic turbulence with periodic boundary conditions is able to produce 
at the same time i) density structures morphologically and statistically consistent 
with the observations; ii) prestellar cores correlated with the gas distribution in 
a way qualitatively similar to the observations and with a value of magnetic field 
strength typically observed; iii) a mass distribution of the same prestellar cores 
that agrees with the observed stellar (and prestellar cores) mass distribution 
over the whole range of stellar masses, from brown dwarfs to massive stars; and 
iv) a star formation efficiency consistent with that in observed molecular clouds 
of only a few percent per large scale dynamical time. 
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6 Conclusions 

The main conclusion from the preceding sections is that the statistics of star 
formation is primarily controlled by supersonic turbulence, rather than by grav- 
ity. Star formation takes place in cold molecular clouds, which are part of a 
turbulent cascade in the interstellar medium. The ultimate energy input to the 
cascade comes from supernovae, with a possibly significant contribution from 
local variations of the galactic rotation curve (density waves). The clouds owe 
their existence to random convergence of the interstellar medium velocity field, 
which creates local density enhancements over a range of scales. The internal, 
supersonic and super- Alfvenic velocity field in molecular clouds is responsible 
for their fragmentation, thus preventing global collapse but triggering local col- 
lapse at the many local density maxima whose mass exceeds the local Jeans' 
mass. Such prestellar cores are formed as sheet corrugations and filamentary 
density enhancements, and are taken over by self-gravity only after they have 
been shaped by the turbulence. 

The velocity field of the cascade is dominated by power in solenoidal (shear- 
ing) motions, even though it is supersonic and super- Alfvenic. Its spectrum of 
kinetic energy is less steep than its velocity and magnetic field power spectra, 
which explains how conditions can be super- Alfvenic on small (molecular cloud) 
scales, even though there is rough equipartition between magnetic and kinetic 
energy density on large (disk thickness) scales. 

A Salpeter like IMF is the result of the near-self-similar, power law nature 
of turbulence in molecular clouds, in combination with density jump amplitudes 
determined by MHD-shock jump conditions. 

Star formation (at least in our galaxy) bites its own tail; it is driven by 
supernovae and at the same time the birth of massive stars gives rise to new 
supernovae that re-enforce the driving. External sources of turbulence, such as 
kinetic energy from galaxy collisions and merging may be the primary driving 
agent in star-burst galaxies. 

Different physical conditions (primarily higher temperatures and lower metal 
abundances) in the Early Universe would lead to higher mass at the low-mass 
cut-off, and a much weaker magnetic field would lead to a steeper IMF slope. 
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